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I. INTRODUCTION 

It is well known that the small- a; behavior of struc- 
ture functions of deep inelastic scattering is determined 
by the hard pomeron contribution. In the leading order 
the pomeron intercept is determined by the BFKL equa- 
tion [T] and the pomeron residue (the 7*7*-pomeron ver- 
tex) is given by the so-called impact factor. To find the 
small- a; structure functions in the next-to-leading order, 
one needs to know both the pomeron intercept and the 
impact factor. The NLO pomeron intercept was found 
many years ago [5] but the analytic expression for the 
NLO impact factor is obtained for the first time in the 
present paper. 

We calculate the NLO impact factor using the high- 
energy operator expansion of T-product of two vector 
currents in Wilson lines (see e.g the reviews [SlH]). Let us 
recall the general logic of an operator expansion. In order 
to find a certain asymptotical behavior of an amplitude 
by OPE one should 

• Identify the relevant operators and factorize an am- 
plitude into a product of coefficient functions and 
matrix elements of these operators 

• Find the evolution equations of the operators with 
respect to the factorization scale 

• Solve these evolution equations 

• Convolute the solution with the initial conditions 
for the evolution and get the amplitude. 

Since we are interested in the small- a; asymptotics of DIS 
it is natural to factorize in rapidity: we introduce the 
rapidity divide r\ which separates the "fast" gluons from 
the "slow" ones. 
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As a first step, we integrate over gluons with rapidities 
Y > r\ and leave the integration over Y < r\ for later time, 
see Fig. 2. It is convenient to use the background field 
formalism: we integrate over gluons with a > a = e' 1 
and leave gluons with a < a as a background field, to 
be integrated over later. Since the rapidities of the back- 
ground gluons are very different from the rapidities of 
gluons in our Feynman diagrams, the background field 
can be taken in the form of a shock wave due to the 
Lorentz contraction. To derive the expression of a quark 
(or gluon) propagator in this shock-wave background we 
represent the propagator as a path integral over vari- 
ous trajectories, each of them weighed with the gauge 
factor Pexp(i<7 J dx^A^) ordered along the propagation 
path. Now, since the shock wave is very thin, quarks 
(or gluons) do not have time to deviate in transverse di- 
rection so their trajectory inside the shock wave can be 
approximated by a segment of the straight line. More- 
over, since there is no external field outside the shock 
wave, the integral over the segment of straight line can 
be formally extended to ±oo limits yielding the Wilson- 
line gauge factor 



U% = Pexp 



Al(x) = d 4 k9(e^ -\a k \)e^A^(k) 



(1) 



where the Sudakov variable a k is defined as usual, k = 
&kPi + fikPi + k±. (We define the light-like vectors p\ 

2 

and P2 such that q = p\ — xbPi and p^ = P2 + r ^ L Pi 
where pn is the nucleon momentum). The structure of 
the propagator in a shock-wave background looks as fol- 
lows (see Fig. 1): 

[Free propagation from initial point x to the point of in- 
tersection with the shock wave z\ 

x [interaction with the shock wave described by the 
Wilson-line operator U x ] 

x [Free propagation from point of interaction z to the 



2 



final point y] . 




FIG. 1. Propagator in a shock- wave background 

The explicit form of quark propagator in a shock- wave 
background can be taken from Rcf. [5 



x, >0>y, 



(2) 



2n 2 (x- z) 4 



■2ir 2 {x- 



where we label operators by hats as usual. Hereafter use 
the notations x + 



[A 



V s + 



(and our metric is (1,-1,-1,-1)). Note that the Regge limit 
in the coordinate space corresponds to x* — >• oo, — > 
— oo while x_L,yj_ are fixed, see the discussion in Refs. 

nam. 

As we mentioned above, the result of the integration 
over the rapidities Y > 77 gives the impact factor - the am- 
plitude of the transition of virtual photon in two- Wilson- 
lines operators (sometimes called "color dipole"). The 
LO impact factor is a product of two propagators ([2]), 
see Fig. [2] 




FIG. 2. Impact factor in the leading order. Solid lines repre- 
sent quarks. 
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x ^^[2(-Ci)(-C 2 )-- 2 (Ci-C2)] + 0(a s ) 
Here we introduced the conformal vectors [El |9] 

yfe ,Pi 2 . v Vs/Pi 2 . \ 
( x p 2 + xj_) - — ( y p 2 + yj_) 



2.x 

Ci = (y + z f±P2 + Z i± ), 



2y* 



(4) 



and the notation 1Z = „, K i^ 2 l % . The above equation 

2(k<i)(k<2) _ ^ 

is explicitly Mobius invariant. In addition, it is easy to 
check that ^|-(r.h.s)=0. 

Our goal is the NLO contribution to the r.h.s. of Eq. 
([7]), but first let us briefly discuss the three remaining 



steps of the high-energy OPE. The evolution equation 
for color dipoles has the form [10] 



2tt 2 



MKuln (5) 



*13^23 

N c tr{U^U^}} + NLO contribution 



(To save space, hereafter Zi stand for Zi± so z\ 2 = z\is_ 
etc.) The explicit form of the NLO contributions can 
be found in Refs. [3J [TTJ Q2] while the agrument of the 
coupling constant in the above equation (following from 
the NLO calculations) is discussed in Refs. ( |13l 114] ). 

The next two steps, solution of the evolution equation 
([5]) with appropriate initial conditions and the eventual 
comparison with experimental DIS data are discussed in 
many papers (see e.g. [15]). 



II. CALCULATION OF THE NLO IMPACT 
FACTOR 



Now we would like to repeat the same steps of operator 
expansion at the NLO accuracy. A general form of the 
expansion of T-product of the electromagnetic currents 
in color dipoles looks as follows: 

(x - y) 4 T{^(^)7^^)^(2/)7>(2/)} 



d 2 Z\ d 2 z 2 



-12 



h h °( 



^[1 + ^tr^^} 



+ jd 2 z 3 I^°(z 1 ,z 2 ,z 3 ;r l ) 

x [tr{t% 0#}tr{C% U%] - N c ti{K U%}] } (6) 

The structure of the NLO contribution is clear from the 
topology of diagrams in the shock-wave background, see 
Fig. [3 below. Also, the term ~ 1 + can be restored 
from the requirement that at U = 1 (no shock wave) one 
should get the perturbative series for the polarization 
operator 1 + ^ + 0(a 2 ). 
In our notations 



n 2 



t 6 (k-Ci)(k-C 2 ) 



dx^dy 



-[( K -Ci)(«-C 2 )-^ 2 (Ci-C2)]. 



(7) 



which corresponds to the well-known expression for the 
LO impact factor in the momentum space. 

The NLO impact factor is given by the diagrams shown 
in Fig. [3j The calculation of these diagrams is similar 
to the calculation of the NLO impact factor for scalar 
currents in N = 4 SYM carried out in our previous pa- 
per [12] . The gluon propagator in the shock- wave back- 
ground at > > in the light-like gauge p 2 A^ = 



3 





FIG. 3. Impact factor in the next-to-leading order. 



is given by [THl E] 



x*>0>y. 



d 4 z 6(z t ) 



(8) 



x x*9-U-P2»{x-z)f uab 1 yj^ - p 2 „(y - z){ 



ir 2 [(x — z) 2 + ie] 



where -i- can be either a 1 . 



which leads to 



the same result. (This is obvious for the leading order 
and correct in NLO after subtraction of the leading-order 
contribution, see Eq. (15) below). 

The diagrams in Fig. |3Ji,b can be calculated using the 
conformal integral 



x 2 y 2 (x — y) 2 



- (i «-> v (9) 



which gives the 3-point iprpF^ Green function in the 
leading order in g. Using Eqs. ([2]), Q and perform- 
ing integrals over z,'s and taking traces one gets after 
some algebra the NLO contribution of diagrams in Fig. 
[3] in the form 



~i Z l, Z 2 ,Z 3 ) = If (Z 1 ,Z 2 ,Z 3 ) +1% (z l7 Z 2 ,Z 3 ) (10) 

where 

ir(zi,z 2 ,z 3 ) 

as rw f ^ ^ (Ci • C2) f°°daiciaz 3 



(11) 



4^ v(zi ' Z2) (Ci-c 3 )(Ci-C3)y oc 



and 

0*2)^(21,22,23) 
„ IV Ca) d 2 



n 2 



16tH» («-Ci)(«-Ca) 



(12) 



(« • Cs) daPdy" L (Ci • Cs) 



( K -Ci) 2 , (K-Ci)(«-Ca) 



+ 



(Ca • Cs) 



(«-Ci)(«-C3)(Ci-Ca) « 2 (Ci • C 2 )l + («-Ca) 



(Ci • C 3 )(Ca • Cs) (Ca-Cs) J ' («-Cs) 2 

a 2 r(«-Ci)(«-Cs) « a (Ci-Cs)- 



c^%* L (C 2 • Cs) 



2(Ca • Cs) 



+ (Ci ^ Ca) 



(recall that z 2 



- 2(G-0) and = ^(«-G))- We 
obtained this expression at a;, > > but from the 
conformal structure of the result it is clear that this ex- 
pression holds true at a;* < < as well. 



The integral over a in the r.h.s. of Eq. ( 11 ) diverges 



This divergence reflects the fact that the contributions 
of the diagrams in Fig. [3] is not exactly the NLO impact 
factor since we must subtract the matrix element of the 
leading-order contribution. Indeed, the NLO impact fac- 
tor is a coefficient function defined according toEq. ^. 
To find the NLO impact factor, we consider the opera- 
tor equation (J6| in the shock-wave background (in the 
leading order (U Z3 ) A = U Z3 ): 

{T{^(x)rMx)kv)i u i>{y)})A 



d 2 z 1 d 2 z 2 LO 



I^(x,y;z 1 ,z 2 )(tT{UlU^}) J 



'12 



d 2 Z\d 2 Z 2 2 7 NLO/ \ 

i d z 3 I {x,y;z l ,Z2,z 3 ;r)) 
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MU Z1 Ul }tv{U Z3 C/t } - N c tv{U Zl Ul }] 



(13) 



The NLO matrix element (T{ip(x)^ip(x)ip{y)'y"i'(y)}) A 
is given by Eq. ( 10 1 while the subtracted term is 



2tt 2 



d 2 z-Ld 2 z 2 LO 



-12 



da 



a 



d z 3 



13^23 



x [tv{U Zl Ul 3 }tT{U Zs Ul} - NMU Zl Ut 2 }} (14) 

as follows from Eq. ([5]). The a integration is cut from 
above by a — e 11 in accordance with th e de finition of 
operators U 71 , see Eq. |l]). Subtracting (14 1 from Eq. 
( 10 1 we get 



V (21, 22, z 3 ; 77) = If "(21,22, 23; 77) + I% v (zi,z 2 ,z 3 ), 
I? v (x,y;zi,Z2,z 3 ;ri) 



" s /-LO/. . \ ^12 

2 n 2 V" \ 2 > J2 ^2 



13^23 



II 



■da 



a 



da 
(i a 



'12 



(15) 



Note that one should expect the NLO impact factor to 
be conformally invariant since it is determined by tree 
diagrams in Fig. [3] However, as discussed in Refs. [4j [7J 
[TT] . formally the light-like Wilson lines are conformally 
(Mobius) invariant but the longitudinal cutoff a < a in 
Eq. (1) violates this property so the term ~ In aZ 3 in the 
r.h.s. of Eq. (15) is not invariant. As was demonstrated 



in these papers, one can define a composite operator in 
the form 



Mu Zl ut 2 }] a 



(16) 
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where a is an arbitrary constant. It is convenient 
to choose the rapidity-dependent constant a — > ae -2 '' 
so that the [trjC/^ &t 2 }] C ° n does not depend on 
i] = In a and all the rapidity dependence is encoded 
into a-dcpcndcncc. Indeed, it is easy to see that 

iMU 2l Ul 2 }] C a 0nf = and £[tr{U Zl Ul 2 }] C a ° nf is de- 
termined by the NLO BK kernel which is a sum of the 
conformal part and the running-coupling part with our 
0(a 2 s ) accuracygKH]. 



Rewritten in terms of composite dipoles (16), the op- 
erator expansion Q takes the form: 

T{^{x)^ijj{x)^{y)^ u ^{y)} 

{Jj?(*l,*)[l + -] MU*M 2 }]a 
z \2 1 " 

+ JePza I^°(z u z 2 ,z 3 ;a) 

x mu z M^{u z M 2 } - NMU Zl Ul 2 }]a} (17) 

We need to choose the "new rapidity cutoff" a in such 
a way that all the energy dependence is included in the 
matrix element (s) of Wilson-line operators so the impact 
factor should not depend on energy. A suitable choice of 
a is given by ao = — k~ 2 + it = — 4z»t>» gQ wg ^^ am 

syx y) 



(x - y) 4 T{^(xh^(x)^(y)Y^y)} 



(18) 
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Here the composite dipole [tr{U Zl Ul 2 }] ao is given by 



Eq. (16) with a = - ^fqfc + *e while 7^q(zi,z 2 ) and 
-^2 "0*i> z 2, 23) are given by Eqs. Q and (12), respec- 
tively. 



III. NLO IMPACT FACTOR FOR THE BFKL 
POMERON 

For the studies of DIS with the linear NLO BFKL 
equation (up to two-gluon accuracy) we need the lin- 
earized version of Eq. (18). If we define 



U a ( Zl ,Z2)=l-^r[tT{U Zl UL}] a (19) 



and consider the linearization 



^tr{U Zl Ul 3 }tr{U z Mj ^r{U Zl Ui 2 }} ao 



U(zi,Z2) -U(zi,za) -~U(z 2 ,z 3 ) 



one of the integrals over in the r.h.s. of Eq. ( 18 ) can 
be performed. The result is 

1 (x - y ) 4 T{^(xh^(x)kyh u Hy)} (20) 



N, 



dn a dn^ f dz\dz- 



dx^ 1 dy v 
where 

^Lo^y; !'^) = ft 



7T 



-7-NLOl 



2 g a/ j (Ci-C 2 )-Ci a C2 3 -C 2 a Cf 
7r«(«-Ci)(K-Ca) 



(21) 



(see Eq. Q and 
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2 
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(22) 

Here one easily recognizes five conformal tensor struc- 
tures discussed in Ref. [T5] , 
While it is easy to see that 

-, — 7 rj -7: — -7— Xx9(x,y;zi) = (23) 

dx^ {x-y) A dx» dy v a(3K ,y ' ; v ; 

one should be careful when checking the electromagnetic 
gauge invariance in the next-to-leading order. The reason 
is that the composite dipole U a °(zi, z 2 ) depends on x via 
the rapidity cutoff ao = 
get 

d 1 dn a 8kP f dz\dz 2 



0^ so from Eq. p) we 



dx^ (x — y) 4 dx^ dy v J z 



U a0 {zx,Z2)X^ o {x,y;zi) 



' {x-yf 



P2 



{x-vY 



dn a Ok 13 f dzxdz % LO d 

—± — x a p \ x ,y,Zi) a—U a (z 1 ,z 2 ) 



dx^ dy v 



(24) 



Using the leading-order BFKL equation in the dipole 
form (linearization of Eq. ([5])) 

d 



-P- 
-12 



a s N c f 2 z\ 
a—U a (zi,z 2 ) = —r^- dz 3 ^- 2 
da 4tt 2 J z( 3 z 2 3 

x [U a (zi,z 3 ) +U a (z 2 ,z 3 ) -U a (zi,z 2 )] 



(25) 
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we obtain the following consequence of gauge invariance 
d 1 



should be the Fourier transformation of Eq. ( 22 1 which 



dXft (x — y) 4 dx^ dy 



v -t- a p yx,y,z z ) 



(26) 



7T 8 x* (x — y) 6 
AnTZ 5 1 \ d 



1 , d In k 2 

211 ' dy v 



[ln( K -Ci)+ln( K -C 2 )] 



We have verified that the expression (22 1 satisfies the 
above equation. 



IV. CONCLUSIONS AND OUTLOOK 

We have calculated the NLO impact factor for the vir- 
tual photons both in the non-linear form ( 18 ) and with 



the linear (two-gluon) accuracy ( 20 ) . Our results are ob- 
tained in the coordinate representation so the next step 



would give the momentum-space impact factor conve- 
nient for phenomenological applications (and available at 
present only as a combination of numerical and analytical 
expressions [TT)]L The study is in progress. 
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